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Properties of bosonic atoms in a small system with a periodic quasi one-dimensional circular 
toroidal lattice potential subjected to rotation are examined by performing exact diagonalization 
in a truncated many body space. The expansion of the many-body Hamiltonian is considered in 
terms of the first band Bloch single particle functions. No assumption regarding restriction to 
nearest-neighbor hopping (tight-binding approximation) is involved. Plots which correspond to the 
zero-temperature phase diagrams of Fisher et al. are obtained and their salient features, in 
remarkable quantitative correspondence with the results available for much larger systems, are dis- 
cussed. The parameter domain associated with the 'Mott insulator to superfluid' transition in the 
ground state corresponds to progressive hopping induced boson delocalization. This is studied in 
terms of the spectral representation of one-body reduced density matrices and thereby linked to es- 
sentially complete dominance of a single eigenstate. Ground state properties relating to superfluidity 
are examined in the context of two-fluid phenomenology. The basic tool, consisting of the intrinsic 
inertia associated with small rotation angular velocities in the lab frame, is used to obtain ground 
state 'superfluid fractions' numerically. They are analytically associated with one-body, uniform 
solenoidal currents in the case of the adopted geometry. These currents are in general incoherent 
superpositions of contributions from each eigenstate of the associated reduced one-body densities, 
with the corresponding occupation numbers as weights. Full coherence occurs therefore only when 
only one eigenstate is occupied by all bosons. The obtained numerical values for the superfluid 
fractions remain small throughout the parameter region corresponding to the 'Mott insulator to 
superfluid' transition, and saturate at unity only as the lattice is completely smoothed out. 



I. INTRODUCTION 

It is a widespread idea that superfluid behavior is a 
consequence of Bose-Einstein condensation. In fact, it 
was shortly after its first observation in experiments with 
Helium that London conjectured this connection, in what 
became known as London's hypothesis In that con- 
text, superfluidity was initially characterized in terms of 
a set of emergent experimental phenomena related to the 
flux properties of the phase II of liquid Helium j2[ , which 
involve, for example, flux without friction through very 
narrow capillaries. Here, Bose-Einstein condensation is 
not of easy experimental characterization, and therefore 
the construction of a scenario capable of revealing more 
about the connections between the two phenomena is 
truly compromised. In addition to experimental difficul- 
ties, the task of approaching these concepts theoretically 
presents forbidding complications in liquid Helium sys- 
tems due to strong interactions which hampers the use 
of known approximation techniques. 

An alternate, simpler, scenario for posing analogous 
questions became established in connection with the 
physics of weakly interacting, cold, trapped bosonic 
gases. Bose-Einstein condensation in such dilute systems 
is well understood theoretically 0^ ■ The refinement of 
experimental techniques, following the pioneering work in 
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achieving the Bose-Einstein condensation in 1995 0, Bj 
allowed for extensive control and manipulation of such 
systems and their use to investigate many properties of 
condensation in both mean-field and strongly correlated 
regimes. Studies of the latter, in fact first considered by 
Fisher ct. al. in 1989 in a different context, in terms 
of the Bose- Hubbard model from statistical analysis 0, 
gained big impulse with the experimental realizations 
of the optical lattices in 2002 [l3] and constitute nowa- 
days a very active part of the research in systems of cold 
atoms [lll,[ll. 

Superfluid behavior related to flux properties in weakly 
interacting Bo se g ases, on the other hand, is not so well 
characterized |13l |. From the theoretical point of view, 
the superfluidity in the Gross-Pitaevskii limit has been 
established rather recently on a basis of a criteria based 
on inertia and two-fluid model arguments [ij]- Much 
of the ongoing experimental research involve setups that 
are capable of investigating superfluid behavior in terms 
of this same criteria, and in particular, to account for ef- 
fects of rotation of containers [H, [Tsl - flTj . This illustrates 
that despite the achievements obtained so far, there is 
still a large class of phenomena related to superfluidity 
in Bose condensates of atomic gases that requires more 
careful analysis. Nevertheless, computational advances 
combined with realizations of sufficient realistic models 
have been allowing for the elucidation of many properties 
related to superfluidity in these systems, like for instance 
the Josephson oscillations [l8| . 

In this work we study properties relating to Bose- 
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Einstein condensation, in the sense of 'macroscopic' oc- 
cupation of a single one-body wavefunction, and to su- 
perfluidity in a weakly interacting Bose gas trapped in 
a rotating quasi-lD ring-shaped Kronig-Penney trap at 
low mean occupation per site, which constitutes a sim- 
plification of the optical lattices implemented experimen- 
tally. Despite being a small system, to which, in par- 
ticular, the concept of a phase-transition strictly speak- 
ing does not apply, it still deploys 'precursor' features 
which can be related even qualitatively to those which 
have been both observed in real systems and supported 
by approximate computational results obtained for con- 
siderably larger systems. In order to study superfluidity 
in terms of the criteria used in , we consider an exter- 
nally imposed rotation ('cranking') of the lattice. From 
this a 'cranked' extension of the Bose-Hubbard Hamil- 
tonian results, which is conveniently diagonalized in a 
single-particle basis of Bloch functions and single-particle 
energies, by calculating explicitly the required two-body 
matrix elements. In the calculations reported here this 
basis is truncated to the first band only. Within the 
bounds set by this limitation, the treatment effectively 
relaxes the constraints inherent to the tight-binding ap- 
proximation which underlies the standard usage of the 
Bose-Hubbard model. In terms of the results of the nu- 
merical many-body diagonalization one can then obtain 
'superfluid fractions' (in the sense of the two-fluid model) 
and associated currents independently of any assump- 
tions regarding a 'condensed mode'. 

Properties relating to condensation are analyzed in 
terms of the establishment of occupation dominance of 
one of the eigcnstates of the one-body reduced density. 
From numerical diagonalization of the many-body Hamil- 
tonian, numerical results for both dominant occupation 
and superfluid fraction are obtained and compared. A 
finite version of the zero temperature phase diagram Q 
is also computed, and it is seen that the closure of the 
Mott lobes happens for the same range of values at which 
occupation dominance essentially attains its peak value. 

The paper is organized as follows: in Sec.|ll]we present 
the cranked version of the Bose-Hubbard model, whose 
ground-state properties are analyzed in Sec. Illll in terms 
of the zero-temperature phase diagrams. In sections IIVI 
and |V] we discuss the properties related, respectively, 
with condensation and superfluidity, and in Sec. IVII we 
present our conclusions. 



II. THE CRANKED QUASI-MOMENTUM 
HAMILTONIAN 

A. Bloch function representation 

The many-body dynamics of cold bosonic atoms in 
external lattice potentials is strongly dominated by the 
two basic ingredients consisting of hopping and short 
range (repulsive) two-body interaction effects. In a tight- 
binding regime hopping is dominated by nearest neighbor 



processes and two-body effects are dominated by 'on-site' 
contributions only. These ingredients are combined in the 
Bose-Hubbard model Hamiltonian 

Hbh = - J ^ (afflj + a]ai) + 

<i,j> 

+ ^^0'\a\aia„ J,U > 0, (1) 

i 

where the first sum is restricted to nearest neighbors 
< i,j > and the operator (aj) destroys (creates) a 
bosonic particle in site i. In order to connect this Hamil- 
tonian involving 'sites' to a more basic description in- 
volving spacial coordinates one associates the sites to 
amplitudes defined in terms of the first band Wannier 
function, which is sufficiently well localized in space in 
the tight-binding regime. As is now well known, competi- 
tion between localization, favored by two-body repulsion, 
and delocalization, favored by the hopping term leads to 
a quantum phase transition between a 'Mott insulator' 
phase and a 'supcrfiuid' phase in the ground state as the 
relative importance of the two parameters J and U of the 
model is varied P, • 

Specializing now to the case of a finite, one-dimensional 
potential array consisting (for convenience) of an odd 
number of sites M with periodical boundary conditions 
(i.e., a ring-shaped one dimensional array of sites), inter- 
esting symmetries become manifest by changing to the 
representation which diagonalizes the hopping term of 
the Hamiltonian ([1]). This is achieved by a (discrete) 
Fourier transform 

1 *^ 

^ n—l 

where q = — , . . . , 0, . . . , ^^-j^. In terms of the new 
'quasimomentum' bosonic operators Aq, the Hamil- 
tonian reads 

Hbh = -2J J] cos(^)iti, 

9=0 

+ ^ E '^A^'^^i + 92 - 93 - qA)A\^Al^Aq^Aq^. 
ij 

(3) 

The Sm{q) in the two-body term is the modular Kro- 
necker delta, equal to one if q is an integer multiple of M 
and zero otherwise. It indicates 'modular' conservation of 
total quasimomentum (Qt = mod (^^qng, M), with 
Uq the number of particles occupying the state of quasi- 
momentum q), i.c, with allowance for Umklapp processes. 
This in fact reduces the many-body Hamiltonian matrix 
to block-diagonal form, each block being associated to a 
value of Qt, which assume the same values of q. The 
presence of tight-binding assumption in this representa- 
tion is manifest in the cosine law for the single-particle 
energies and in the single common value U/M for the 
two-body matrix elements. 
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Limitations of the tight binding regime which are buih 
into the Bose-Hubbard Hamiltonian can be Ufted without 
substantiahy increasing eventual computational costs by 
first returning to a configuration space representation of 
the Hamiltonian. We do this by keeping the overall ge- 
ometry consisting of a regular toroidal array with mean 
radius R and containing M angular domains separated 
by potential barriers, with a one-dimensional lattice con- 
stant given by Ic = 2-kR/M. Furthermore, we assume an 
effective one-dimensional regime in which the transverse 
amplitude is independent of angle and effectively frozen 
in its ground state. In this context the relevant Hamilto- 
nian can be written in terms of angle dependent bosonic 
field operators i)[(p)^ i^Hf) as 



H 



dip V'^(ly3) 



A 



Vlattiv) 



(4) 



in which m is the boson mass and Viatti^) is the exter- 
nal lattice potential, accounting both for tight transverse 
confinement and for the angular periodicity. The effec- 
tive one-dimensional strength parameter A is related to 
the strength of the usual effective two-body contact in- 
teraction given in terms of the s-wave scattering length 
a, A = ATrti^a/m, as 



A = A 



\Mr±)\' 



Here u'o(r^) is the frozen transverse amplitude in the 
array. 

Before making contact with the Bose-Hubbard form 
(U) of the Hamiltonian, we consider a further extension 
of the effective one-dimensional Hamiltonian (|4]) to in- 
clude uniform rotation of the lattice with angular veloc- 
ity uj around the axis of the toroidal structure. By trans- 
forming to the reference frame rotating with the lattice 
potential, the required effective Hamiltonian becomes 



dip ip^ (ip) 



1 fh d 



\i dip 



A 



muR^ 



(5) 



1p{ip). 



This differs from Eq. (j4]) just by the replacement 



h d 

i dip 



{I2 — muiR) 



h d 

i dip 



mujR. 



(6) 



One can now make contact with the Bose-Hubbard 
form ([3]) by representing H^^ in the (truncated) single 
particle basis consisting of the M (first band) Bloch 

functions {^^"^((ys)} which diagonalize its one-body part. 
Note that these functions are labeled by the quasimo- 
mcntum q and depend on the cranking angular velocity 



w. This is done by expressing in terms of the bosonic 
operators 



A. 



(^) 



A. 



which leads to 



d^^<^\^)^{^) and 
d^ ((^)^,((^), 



M-l 



J2 e(-)i(")U(- 



q=0 

A 



(7) 



V 4('^)t4('^)t4('^)4(") 

A -rlq-^ -ri-nn -^o-i -^n 



^92 



^93 ^94 



Here the e^'' are the (first band) single particle Bloch 

eigenvalues and the objects ^^j^'j = 2^g".g2.93,94 ^'^^ inte- 
grals over products of four Bloch functions which, to- 
gether with the strength parameter A, constitute the re- 
quired two-body matrix elements in the adopted single- 
particle basis. The single particle energies are now the 
Bloch eigenvalues which, as shown in continuation, have 
spacings consistent with the cosine rule of ([3]) in the tight 
binding regime; the two body matrix elements may and 
indeed do increasingly show fluctuations as one leaves the 
extreme tight binding limit. 

Finally, it is worth mentioning explicitly that trunca- 
tion of the single-particle basis to the first band Bloch 
states, together with the modular conservation of total 
quasi- (angular)momcntum in ([T]) implies that there are 
no 'mean field' two-body effects, i.e. all two-body ma- 
trix elements involving just one-particle changes of state 
vanish. As a consequence, reduced one-body density ma- 
trices of non-degenerate stationary many-body states of 
([7]) will be diagonal in this basis. 



B. Numerical features 

In order to implement this many-body model numer- 
ically we first adopt a schematic realization of the lat- 
tice potential in terms of the Kronig- Penney model with 
(5— function barriers 



Viatt -^i^5{ip~ —{n 

Tl = 



1)), 



where 7 is the strength parameter for the (5— function 
barriers. For this choice, Bloch functions are given as 
analytic expressions involving few numerical parameters 
which are easily obtained, together with the Bloch eigen- 
values, by solving transcendent algebraic equations nu- 
merically. 

The behavior of the Bloch eigenvalues as a function 
of the cranking angular velocity w, given in units of 
Wo — h/mR^ (energy being measured in units of eo = 
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FIG. 1. First band Bloch single-particle energies (in units of 
eo = h^/mR^), for M = 5, 7 = leo (a) and lOOeo (b), as 
functions of the cranking angular velocity lo (in units of cuo; 
(c): Bose-Hubbard single-particle energies per unit hopping 
parameter also as a function of to. Note the relative com- 
pression of lower levels in (a). One quarter of the width of 
the Bloch band may be used as an effective Bose-Hubbard 
hopping parameter J for each value of 7, as discussed in the 
text. 



h^/mR^), is shown for two values of the barrier param- 
eter 7 in Fig. [T] (a) and (b). Part (c) of the same figure 
shows the uj dependence of the tight-binding regime one- 
body eigenvalues as implemented in the Bose-Hubbard 
model, given by the function — 2cos(|j((j' — uj)) (cf. Eq. 
dS]). While this cosine function is clearly capable of repro- 
ducing the relative spacings of the Bloch single-particle 
energies very accurately for sufficiently large values of 7, 
there is a clear low-energy bunching of the single particle 
levels for the lower value of 7. In the following we use 
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FIG. 2. Typical behavior of the two-body matrix ele- 
ments. The notation (91,52,53,54) denotes the integral 
J dLp(l}*^{tp)(t>*,^{(p)(j>g^{ifi)(j>^^{if). Unlike strictly quasi-mo- 
mentum conserving matrix elements, Umklapp matrix ele- 
ments vanish together with barrier strength parameter at 
J = 25/32eo. 



this to define a parametrization of the barrier strength 
in terms of an effective hopping coefficient J as 



J 



(8) 



where and are respectively the upper and the 
lower bounds for the first band cranked Bloch eigenvalues 
{e^"^} for a given value of the barrier strength parameter 
7. For simplicity we refer to the numerator of Eq. ([8]) as 
the 'energy width of the first band'. This coincides with 
the Bose-Hubbard definition of the hopping parameter in 
Eqs. Q and ([3|) in the tight-binding domain. Note also 
that the dependence of the Bloch energies with the crank- 
ing angular velocity leads to single-particle energy level 
crossings at integer and and half-integer values of uj/uiq. 
Since the quasimomentum is a good quantum number for 
the individual Bloch states, these level crossings will have 
a decisive role in determining the total quasimomentum 
of cranked many-body ground state, as discussed in the 
following section. 

The four wavefunction integrals Iqijq2,q3,qi are easily 
evaluated in terms of the Bloch functions. Some sample 
results are shown in Fig.[2J for M = 5 and a; = 0, as func- 
tions of the effective hopping parameter J just defined. 
Salient features here are the complete bunching on a sin- 
gle value To in the tight-binding limit J ^ {Iq ^ 0.239 
in the present case), and the strong quenching of Umk- 
lapp matrix elements in the opposite limit of very large 
hopping, as compared to angular momentum conserving 
matrix elements. Note that the definition of the effec- 
tive hopping parameter leads to J/eo — > 25/32 = 0.78125 
as the barrier strength parameter 7 — > 0, for the present 
case with M = 5. In order to facilitate comparison with 
results obtained using the Bose-Hubbard Hamiltonian 
parameters as in ([1]) and ([3]), we again use the tight- 
binding limit to define an effective two-body parameter 
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U to replace the two-body constant A as 

U = IqMA (9) 

in terms of which the two-body part of eq. (O coincides 
in the tight-binding hmit to that of the Bose-Hubbard 
model, eq. (|3]). 

III. THE GROUND-STATE PHASE DIAGRAM 
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FIG. 3. Particle addition energies as functions of J for uj = 
and to — 0.48 using (the) single-particle energies, and with 
the effective hopping parameter defined in the text. Here 
Uo — 0.06eo, M — 5 and A'' varies from to 11 (from the first 
to the last line going upwards). In (a) the closure of the Mott 
lobes is seen to happen in the neighborhood of J/Uo = 0.2, 
while in (b) it is affected by the angular velocity and happens 
for higher values, around J/Uq — 0.6. 

The physical properties of systems described by the 
Bose-Hubbard Hamiltonian ([3]), in particular, the dif- 
ferent dynamical regimes prevailing in different param- 
eter domains have been studied extensively using di- 
verse approximation schemes and/or computational tech- 
niques. The focus of these studies is the thermodynamic 
limit in which the number of sites and the number of 
bosons go to infinity, at fixed finite mean occupation per 
site. The one-dimensional case has been treated rather 
recently by Kiihner, White and Monien [l^, including, in 
particular, the 'ground state phase diagram' Q in which 
different phase domains are identified in a ^/U x J/U 
diagram, fi being the chemical potential. 



Our purpose in this section is to indicate how one can 
obtain phase-diagram information by using the results of 
exact many-body diagonalization of the Hamiltonian ([7|) , 
albeit with feasibly small number of sites and of bosons. 
Under these circumstances there will of course be no case 
for taking thermodynamic limits and determining precise 
phase boundaries, but as will be shown explicitly there 
arc clear 'precursor features' in the small number solu- 
tions which clearly identify, surprisingly even quantita- 
tively, several thermodynamic limit properties. 

Since in this context we deal always with systems hav- 
ing a fixed, given number of particles, a replacement must 
be devised for the chemical potential fi. We thus replace 
fi by the addition energy Ajv, N being the number of 
bosons, defined as A^r ~ Eo{N -t- 1) — Eo{N) — e, where 
Eq (N) denotes the ground state eigenvalue for N bosons 
in the chosen number of sites, and e is the average of the 
single particle (first-band Bloch) energies. Note that e 
vanishes for the cosine law of the Bose-Hubbard Hamil- 
tonian ([3]). With this replacement, the axes of the graph 
corresponding to the ground state phase diagram become 

M Aat _ EojN + 1) - EoiN) ^ e J_ 

U U U ' U ^ ' 

where we use the quantities U and J as defined in Eqs. ^ 
and ([HI) in order to characterize the dynamical parameters 
of the many-body Hamiltonian ([7]). 

Results obtained for il/ = 5 sites and = to 11 
bosons are shown in Fig. [3] (a), for the usual situation of 
vanishing cranking angular velocity, a; = 0, the strength 
of the two-body effective interaction having been fixed 
at A = Q.QQTi' /mE? . The value of the barrier strength 
parameter 7 has been varied to cover the desired range 
of values of J/U . Actually the whole parameter domain 
covered by this graph falls within the tight-binding do- 
main in which the results obtained using the Hamilto- 
nian ([7]) differ very little from what one obtains using 
the Bose-Hubbard Hamiltonian ([3]) itself [l^l • Each curve 
in this graph shows the dependence on J/U of the {U- 
scaled) addition energy A at for one of the values of A^. 
As seen, the bunching of these curves at integer values of 
the JJ-scaled addition energies in the limit J/U — > gives 
room for the Mott Insulator lobes. The shape and range 
of these lobes reproduces even quantitatively the results 
obtained towards the thermodynamic limit, cf. Ref. (loj . 
including in particular the 'reentrant behavior' plj (see 
also Ref. [l^] ) characteristic of the one-dimensional Bose- 
Hubbard model. Moreover, one can show perturbatively, 
in the large J/U limit past the insulator lobes, that the 
spacing between consecutive Aat /U curve s ap proaches 
the value 1/M, M being the number of sites [20|, suggest- 
ing therefore a definite A^-cleavage as one moves toward 
the thermodynamic limit. 

Effects of cranking are illustrated in Fig. [3] (b), which 
differs from (a) in that here a; = 0.48a;o, i.e. just 
shortly before of the first single-particle level crossing 
(see Fig. [T|) . The broadening and lengthening of the Mott 
insulator lobes which is visible in this case is governed es- 
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sentially by the Bloch energies modified by rotation. This 
effect can be understood in terms of the effective reduc- 
tion of the kinetic energy (cf. Eq. ([6])) which quenches 
hopping thus favoring the insulating pliase. 



IV. PROPERTIES RELATED TO 
CONDENSATION 

All the information relevant for one-body observables 
contained in the correlated many-body ground state 



I* 



) of the Hamiltonian ([7]) is conveniently extracted 
by evaluating its reduced one body density defined as 

EE ($[,"'^|V>^((p')^(<^)|$o"^)- This is a hcrmi- 
tian kernel of trace equal to the number N of bosons 



I*, 







hence with a purely discrete spectrum. The 



solutions of the eigenvalue problem 

dip p'^^\(p,ip )Xu{^ ) = n^Xuiv) 



lead to the so called natural orbitals x^i'P)^ the asso- 
ciated ei gen values being the corresponding occupation 
numbers j23| . Here, in particular, the fact that the many- 
body states l^o"^^) have good total quasi-momentum, to- 
gether with the adopted truncation to single boson states 
of the first band only, ensures that the associated reduced 
one-body densities are diagonal in the quasi-momentum 
representation, i.e. 



(11) 



which leads to the spectral decomposition p^^\ip,ip ) = 

'^qiT'qCj)^'' {ip)(t>'l^\(p ), with thc Bloch functions ^g"' 
themselves playing the role of the natural orbitals of the 
system. 

Interest here in these objects stems from their rela- 
tion to the characterization by Penrose and Onsager [23| 
of Bose-Einstein condensation in interacting many boson 
systems. It is done in terms of the macroscopic occu- 
pation of a natural orbital, which then plays the role of 
'condensate wavefunction'. Barring important limiting 
procedures which have to be used when dealing with ex- 
tended systems, and notwithstanding non-existence theo- 
rems for Bose-Einstein condensation as a phase transition 
in one dimensional systems [22| . in the present case of 
fixed, finite number of bosons the behavior of the natural 
orbital occupation numbers with the parameters of the 
model Hamiltonian is still revealing of condensation ori- 
ented one-body properties of the correlated many-boson 
state. 

Typical behavior for occupation numbers of the sys- 
tem's ground state natural orbitals is shown in the plots 
of Fig. m (a), (b) and (c) for a situation with commen- 
surate filling (the number of bosons A'' being an integer 
multiple of the number of sites in the array, M), and 
in Fig. m (d) for an incommensurate case. In all cases 
the two-body interaction parameter has been fixed at 
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FIG. 4. Eigenvalues of the reduced density matrix for Uo = 
Omn'^/mR'^ and M = 5. In (a), (b) and (c) po = 1, while 
in (d), Po — 1.2 . In (a), co — while in (b), co — 0.48 and 
in (c), a; = 0.52, i.e. just before and after the single-particle 
level crossing at w = 0.5. At lj = 0.5 there is a migration of 
population from the state with quasi-momentum g = in (a) 
and (b) to the one with q — 1 in (c); (d) illustrates qualitative 
differences in the limit J — >■ between the incommensurate 
and the commensurate situations. 
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A = 0.06h^ /mR"^ . Since in the incommensurate case 
one may have to deal with different values of total quasi- 
momentum due to the w-dependence of the Bloch ener- 
gies, the two cases are better discussed separately. 



We start with the commensurate situation where, re- 
gardless of the values of uj, the ground state of the many- 
body system is always found in the subspace of total 
quasi- momentum Qt = 0. As shown in Figs. U] (a), (b) 
and (c), evaluated for the case N = M, in the extreme 
tight-binding limit J — > 0, one has a fully degenerate 
situation in which each of the orbitals has unit occupa- 
tion, implying that the reduced one-body density is an 
incoherent superposition with equal weights of one body 
densities constructed from each of the first band Bloch 
functions. This is in fact a consequence of the localization 
of each particle in one of the sites, achieved through prob- 
ability amplitudes involving the (Wannier) coherent su- 
perposition of the Bloch functions. In the commensurate 
case po = 1 the localization causes the two-body energy 
to vanish in the limit J — > 0. Increasing the effective hop- 
ping parameter J by lowering the barrier strength leads 
eventually to a situation in which the trace of the re- 
duced one-body density is carried by essentially a single 
eigenvalue, the associated natural orbital being the Bloch 
state with lowest single-particle energy at the considered 
value of the cranking angular velocity w. This implies 
that the spectral decomposition of the reduced one-body 
density essentially reduces to a single term, which corre- 
sponds to the underlying many-body state approaching 
the simple form of a product state in which all bosons are 
coherently delocalized in the collectively occupied Bloch 
wavefunction. Parts (b) and (c) of Fig. 0] illustrate the 
change of role of Bloch orbitals near the level crossings 
at UJ — O.Swq. 



Fig. m (d) shows the evolution under increasing ef- 
fective hopping parameter J of the eigenvalues of the 
reduced one-body density in the incommensurate case 
TV = 6, M = 5 for = 0. Unlike in the commensurate 
case, here one no longer has occupation number degener- 
acy for the natural orbitals at the extreme tight-binding 
limit J — > 0, but an occupation enhancement of the low- 
est Bloch state which in fact typically exceeds the con- 
tribution of the extranumerary boson, an effect which 
may be traced to the action of bosonic enhancement fac- 
tors [l^l . The delocalization process induced by hopping, 
associated with progressive dominance of a single Bloch 
state, is however maintained. 



Finally, it is worth noting that, in all cases, the scale 
over which the transition to delocalization occurs, in 
terms of J/U coincides with that which may be asso- 
ciated with the extension of the Mott insulator lobes in 
the ground-state phase diagram, sec Fig. |31 suggesting 
the connection of these two features. 



V. PROPERTIES RELATED TO 
SUPERFLUIDITY 

The dependence on w of the energy (in the rotating 
system) of the many-body ground state of the cranked 
Hamiltonian ^ , for small values of the angular velocity, 
allows for the determination of an incrtial parameter I 
through the relation 

(N) = {N) + + 0{uj^), w ^ 

which may be related to the 'superfluid fraction' /g, de- 
fined within the phenomenological framework of the two- 
fiuid model of superfluid behavior, as (cf. Ref. Q) 



-L Ring 



NmR'^' 



(12) 



This bears the understanding that, in stationary state, as 
viewed from the rotating frame, the inertia corresponds 
to that part of the fluid which remains stationary in the 
lab frame (the 'superfluid'), while the 'normal' compo- 
nent is carried around with the externally imposed angu- 
lar velocity. 

The inertial parameter I can be expressed as 



1=2 



= i($('^)|^|$M) 
UJ duj 



(13) 



where |$(")) is the state vector for the A^-body ground 
state at cranking angular velocity uj and where the well 
known Feynman-Hellman relation has been used. This 
expression involves the expectation value of the deriva- 
tive of Hi^, Eq. ([U, with respect to ui, which is in fact a 
one-body (albeit nonlocal) operator: 



dH^ 
duj 



dip^^{ip) ( — mR^uj] '0(iy9). 

^ i dip J 



An interesting expression for the matrix element appear- 
ing on the expression (jl3p for the inertial parameter I 
can thus be obtained in terms of the spectral decom- 
position of the reduced one-body density, p^^\ip^ip ) = 

Y.uT^vXv{'^)xl{'^')- It reads 



1 = 



mR^ 



(14) 



where the current jui'f) is given by 



Im (x*{v) — 
mR"^ \ dip 



To obtain this expression, matrix elements involving a 
(ys-derivative have been computed with the prescription 



dip 



dip / dip' ij'^{ip)6'{ip-^')'4>{ip'). 
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The current j^^ {(p) is therefore in general an incoherent 
sum of currents associated to each natural orbital, full 
coherence resulting only in the limiting case of full occu- 
pancy of a single natural orbital. Moreover, it is easy to 
see that the current Jw(<p), together with the (diagonal 
part of the) one-body reduced density, obeys a conserva- 
tion law which, for the stationary states of the cranked 
Hamiltonian ([5]), in fact makes it independent of if (i.e., 
solenoidal) as a result of the time independence of the 
associated one-body reduced density. 

In order to see this, we start from the expression of the 
many-body Hamiltonian in terms of the field operators, 
Eq. ([7]), and use the stationary character of |v]/('^)) 



|*(t) 



to write the vanishing time derivative of the reduced one- 
body density as 



1 

ih 



1$ 



By computing the cp-derivatives with the prescription 



dip^ 

dip I dip ^\p)5 (ip-p )4){p ), 



the commutator can again be expressed in terms of the 
reduced one-body density expressed in terms of its nat- 
ural orbitals with the result 



d X ^ 



Im 



dip 



= 0, 



(15) 

which states the p independence of the current j^j ■ 

It is worth noting that these formal results can be ex- 
tended in a straightforward way to a more general con- 
text, involving stationary many-boson states in a three- 
dimensional external trap cranked around a fixed axis, in 
which case the inertial parameter is given in terms of the 
volume integral of the appropriately weighted tangential 
component of a solenoidal current density. 

Using (jlSp . the expression for the inertial parameter 
I, Eq. dU]), reduces to 



2 27r 

I = mR — 

U! 



Correspondingly, the 'supcrfluid fraction' fs is expressed 
in terms of this current as 



= 1— ■ 



(16) 
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FIG. 5. Superfluid fraction for the system with N — M = 5 
and A = 0.06h^ /mR^ . The curve corresponds to Eq. H12|) and 
the points are obtained using Eq. (fT6|l . 



In fact, the last two factors clearly amount to the time 
integral of the solenoidal current ji^ for the duration of 
one period of the cranking. 

A feature worth pointing out explicitly of the evalu- 
ation of supcrfluid fractions by means of Eqs. (fT2|) and 
(|16p is that it involves no assumptions concerning the 
existence of a condensed mode or of any particular cur- 
rents associated with it. Note that the current has 
itself been obtained as a mere byproduct of the two-fluid 
definition of the supcrfluid fraction. 



A. Numerical results 

In order to deal with the implied limit w — >■ in the 
numerical evaluation of superfluid fraction values from 
either Eq. (|12p or Eq. ([TS]) , we use the fact that the value 
obtained is very stable against variation of the value of cj 
even by orders of magnitude down from oj ^ Q.lfi/mR?, 
as long as the differences involved in obtaining either 
T or are sufficiently above the limitations set by 
machine precision. In view of Eq. (jl4p . this amounts 
to a numerical verification that the w-dependcncc of 
£;('^)(7vr)-£;(o)(iV) is quadratic, while that of j„ is linear 
within such a range. 

The most striking result of the numerical evaluation 
of the supcrfluid fraction for the model system on hand 
is the 'prima facie' absence of any dramatic features in 
the parameter domain which supports the delocalization 
transition and the closing of the Mott insulator lobes (cf. 
Figs. 4 (a) and 5). In fact, unlike the saturation in 
occupation of the dominant orbital state, the superfluid 
fraction remains rather small in this region, showing a 
smooth monotonically increasing behavior on a deflnitely 
larger scale of J values. Fig.[5]shows the agreement of re- 
sults obtained from Eq. p^ . the inertia X being obtained 
numerically from the lo dependence of the many-body 
eigenvalue, and Eq. p6p . the current ju, being obtained 
from the Bloch functions and from the eigenvalues of the 
reduced one-body density. 
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(a) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

J/U 



(b) 
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10 
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11 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 
J/U 

FIG. 6. Quenching of the superfluid fraction at commensurate 
fiUings, TVi = M in (a), N2 = 2M in (b). In both cases curves 
for Ni — \, Ni and A'^i + 1 are shown, and the lowest curves 
correspond to the commensurate fiUing situation. Compare 
the extension of the quenching interval to the closing of the 
corresponding Mott insulator lobes in Fig. [3] (a). 

Effects related to the presence of the Mott insulator 
lobes in the ground-state phase diagram of Fig. |3] (a) do 
appear, however, in the corresponding intervals of J/U 
and consist of a quenching of the calculated superfluid 
fraction in such intervals at the corresponding commen- 
surate filling, relative to the values obtained for incom- 
mensurate filling. This is shown in Figs. [5] (a) and (b) for 
the cases N/M = 1 and 2 respectively, and may thus to 
be associated with the stronger loss of coherence of the 
current (see Eq. (jl4p ') occurring in the case of com- 
mensurate filling. Note that the calculated values of the 
superfluid fraction are small in the whole J/U domains 
associated with the Mott insulator lobes, so that these 
effects are not conspicuous on the scale used in Fig. [51 



falls within the tight-binding domain in which the Bose- 
Hubbard model is an excellent approximation to the 
present treatment of the Hamiltonian ([7|) , and that the 
potential barriers themselves constitute an important 
mechanism coupling the externally imposed rotation to 
the dynamics of the many body system. From this point 
of view the obtained values for the superfluid fraction 
based on the two-fluid model arise from a current, in 
the rotating frame of reference in which the lattice po- 
tential is at rest, that still manages to flow through the 
hindrance of the potential barriers. 

There are two features of this current that deserve 
some remarks. First, the fact that in the present one- 
dimensional circular geometry it is independent of po- 
sition on the circle (see Eq. (|15p ) is clearly a necessary 
consequence of stationarity, which in particular requires 
a time independent one-body density. Moreover, the re- 
striction to the Bloch orbitals of the flrst band, together 
with the (modular) conservation of quasimomentum by 
the Hamiltonian ([7|), imply the vanishing of mean field 
effects of the two-body force, thus causing the first band 
Bloch states themselves to play the role of eigenstates of 
the reduced one-body density. From the one-body equa- 
tion leading to them it follows that each of the individual 
square brackets in Eq. p5|) vanishes for the numerical 
calculations reported here. It should be kept in mind, 
however, that this is just an artifact of the adopted trun- 
cation of the single particle basis, and not true in the 
general case in which, according to Eq. only the 

appropriately weighted sum of them is independent of ip. 

The second feature deserving comment is the fact that 
the current associated with the superfluid fraction is gen- 
erally given as an incoherent sum of contributions as- 
sociated with each of the populated eigenstates of the 
reduced one-body density. Thus, to the extent that it 
favors the occupation of a single such eigenstate, the de- 
localization transition associated with the closing of the 
Mott insulator lobes effectively promotes the complete 
coherence of the full current. In the context of the calcu- 
lations reported here, the effects of incoherence stem from 
the different values of the ((/j-independent) contributions 
of the first band Bloch states together with the respec- 
tive occupation numbers. The special role played by the 
Bloch states is however also an artifact of the adopted 
truncation of the single-particle basis. In general, fur- 
ther many-body dynamical effects will be encoded in the 
structure of the eigenstates of one-body density. 



VI. CONCLUSIONS 



B. Discussion 

The relatively small values obtained for the super- 
fluid fraction throughout the range of barrier strengths 
in which the delocalization transition takes place is per- 
haps not surprising if one keeps in mind that this range 



We presented a numerical study of properties of a small 
system of identical bosons in one-dimensional circular lat- 
tice undergoing rotation. The calculations were based on 
the diagonalization of the standard many-body Hamilto- 
nian in the external lattice potential with contact two- 
body effective interactions truncating the involved single 
particle basis to the first band Bloch states of the ro- 
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tating lattice. With the appropriate scahng of the sys- 
tem parameters, results agree quantitatively with those 
obtained from the simple Bose-Hubbard model in the 
tight-binding limit, but will elsewhere incorporate addi- 
tional effects (such as modified single-boson energy spac- 
ings and hopping effects other than nearest neighbors), 
subject to the limitations imposed by the adopted trun- 
cation of the single-particle orbitals. 

Results obtained for a reinterpretation of the (grand- 
canonical) ground state phase diagram in terms of 'sep- 
aration energies' obtained by comparing ground state 
energies of systems with different numbers of particles 
reproduce results of known grand-canonical calculations 
quantitatively even for a small lattice, including the 
'reentrant behavior' of the Mott insulator lobes charac- 
teristic of one-dimensional systems. The closing of the 
lobes is clearly associated with the delocalization of the 
bosons over the lattice which results from the progressive 
occupation of a single eigenstatc of the one-body reduced 
density matrix. 



Numerical results for ground state supcrfluid fractions, 
obtained in the context of two-fiuid phenomenology, link 
the closing of the Mott insulator lobes with increased 
coherence of a uniform, solcnoidal single particle current 
associated with the intrinsic inertia of the superfiuid frac- 
tion. This association can be formally extended to three 
dimensional traps, also independently of assuming the 
existence of a condensed mode or any particular currents 
associated with it. 

The value of the superfiuid fraction remains small 
throughout the parameter domain associated with the 
transition, and grows monotonically saturating at the 
value 1 only with the complete smoothing out of the lat- 
tice. This suggests that the 'Mott-insulator to superfiuid 
transition' may be a transition related rather to proper- 
ties of coherence and delocalization, than to supcrfiuidity 
in terms of flow propensities in the system. 

FP thanks Jonas Larson and Jani-Petri Martikainen 
for useful comments and discussions. Work has been 
supported by FAPESP and the Swedish Royal Council 
(Vetenskapenradet) . 
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